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This paper deals with the output improvement of heating and cooling cycle by using the work-fluid including 


phase change material. The experimental study is carried out by heat exchange between work-fluid and heat 


transfer surface. The work-fluid is flown to a high temperature or a low temperature heat transfer surface from the 


narrow path. In order to increase the amount of the heat transmission, a trace of Diethylether (boiling point 


34.8 °C), as a phase change material (PCM), is added to the work-fluid. The parameters of the experiment are ad- 


ditive amount of PCM, the rotational speed of the displacer piston and the temperature of heat transfer surface. It 


is clarified that the increasing of engine cycle output is brought by the PCM addition. The effect of PCM addition 


is evaluated by output ratio which is defined from the experimental cycle output data. The requirements for ac- 


quiring the increasing effect of output by adding PCM are clarified. 


Keywords: Phase change material, Heat transfer, Output improvement, Low temperature difference engine 


Introduction 


Low temperature (less than 100°C) heat as factory 
exhaust heat or hot spring heat is almost all thrown away. 
Most of the low temperature heat energy is not used at 
present. More than 3000 hot spring resorts with accom- 
modation [1] exist in Japan. Although the electric power 
generation by using the binary cycle [2] or the thermoe- 
lectric conversion element [3] is beginning to be tried at 
some places, but the construction of equipment involves 
large costs. The installation of expensive large equipment 
is needed to gain the output from the low temperature 
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heat energy. Enlarging of equipment is brought by ex- 
pansion in the heat exchange part and so on. The long 
generated electric output selling period is needed for cost 
recovery. 

The final purpose of this research is the development 
of low temperature difference drive engine supposing use 
in a hot spring resort as a motor for small electric power 
generation. In order that travelers can enjoy seeing it, the 
number of rotations of the engine is established at low 
speed. An engine cycle of this study is aimed at the Stirl- 
ing cycle which can maintain high efficiency also in 
small size. This kind of engine has simple structure, it 
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Nomenclature 


a Coefficient y Volume (m°, mL) 

Aaa Mass additional ratio to Air (mass%) W Output power (J/cycle) 

b Coefficient Ww Output power ratio (-) 

Cv Specific heat (J/(kg:K)) Greek letters 

D Depth (m) v Kinematic viscosity (m/s) 

G Mass (kg) 0 Rotation angle (rad) 

H Height between displacer piston and heat 1 Thermal conductivity (W/(m-K)) 
transfer surface (m, mm) 

L Length (m) Subscripts 

N Number of revolutions (r/min) 0 Initial 

n Number of times (-) 1 heating 

P Pressure (kPa) 2 cooling 

AP Pressure difference (kPa) air Air 

q Exchange heat per time (J/s) Add PCM including 

Q Exchange heat on heating process (J/cycle) b Boiling point 

S(=W,xD) Area of heat transfer surface (m?) cycle Cycle process 

St Stroke (60 mm) d Diethylether 

t Time (s) DP Displacer piston 

T Temperature (°C) i Data step 

AT Temperature difference of heat source (K) mc Mean 

Temperature difference ratio (-) pr Process 

u Mean velocity on heat transfer surface (m/s) w Work-fluid (Air, Air+PCM) 


brings low cost, and it is easy to perform maintenance. 
However, it is difficult to obtain enough output by this 
type of engine, because of its low temperature difference. 
Authors found out the increasing effect of engine output 
by adding a little phase change material which has the 
boiling point between the source of high temperature and 
the source of low temperature in a work-fluid [4]. This 
paper deals with increasing of engine cycle output by 
using work-fluid including phase change material. The 
work-fluid is flown into the heating or cooling surface 
from the narrow path on cycle process. Some researches 
[5-7] on heat transfer evaluation of a work-fluid have 
been published, but these studies were examined by us- 
ing a single phase gas as work-fluid. In order to increase 
the cycle output, a little amount of Diethylether (C4H100, 
molecular weight 74.12 g/moL, density 0.7134 g/cm’ on 
liquid phase, boiling point 34.8°C, colorless liquid), as a 
phase change material, is added to base work-fluid (air). 

This study indicates the result that estimated stable 
drive cycle output experimentally. The experimental pa- 
rameters are the number of rotations 5-20 r/min, the heat 
source temperature difference 45-60K and the mass 
composition ratio of phase change material 0-4mass%. 
The evaluated result is indicated as useful experimental 
correlation equation. 


Experimental Procedure 


Experimental apparatus 


Fig.l shows the schematic diagram of the experimen- 
tal apparatus of this study. The experimental apparatus 
consists of a test section, an actuator, a data logger and 
the 2 set of constant temperature bathes. The test section 
is made of the clear acrylic resin, and thermal insulation 
panels 10 mm in thickness have been stuck on the inter- 
nal surface. The inner size of the test section is 120 mm 
in height, 300 mm in width and 150 mm in length. The 
upper surface and the bottom surface of the test section 
are made of copper plate 3 mm in thickness, and are a 
cooling area and a heating side, respectively. The cooling 
and heating surface temperature are measured by 6 points 
of K type thermocouples that have 0.18 mm diameter. 
The heating and cooling surface temperature are con- 
trolled within +0.2K of the experimental condition. The 
pressure change in the experiment is measured by the 
pressure sensor (Keyence AP-43, range 0~ +1.000 MPa, 
resolution 0.1 kPa, response time | ms) that is attaching 
to the acrylic test section wall. Putting the phase change 
material into the test section is performed using digital 
burette (Continuous E 2.5 mL/rev., minimum scale 0.01 
mL, accuracy +0.2%). A displacer piston that is set in the 
test section is connected to an actuator (Misumi single 


axis robot RSD2 position accuracy +0.02 mm, maximum 
load 25 kg, stroke 50-300 mm) through a SUS seal rod 
(3 mm in diameter), and is moved up and down by com- 
puter control. The increasing of pressure in the test sec- 


tion is caused by moving of displacer piston to upper side. 


The upper side work-fluid is flowed to bottom side by 
raise of displacer piston. The pressure of the test section 
is increased by heat exchange between work-fluid and 
hot heat transfer surface. The decreasing of pressure is 
brought by moving of displacer piston to bottom side. 
The heat exchange cycle experiment is carried out on 
each condition. Work-fluid is passed through 3 mm gaps 
(length of 150 mm) of side wall by up-and-down motion 
of a displacer piston. The process distance of a displacer 
piston is 60 mm. 


Thermophysical properties of test work-fluid 


The work-fluid of this study is air or air-diethylether 
mixture. Diethylether is selected as the phase change 
material which has the boiling point in temperature range 
between heating and cooling surface temperature. Table 1 
shows the thermophysical properties of air and diethy- 
lether. For example, the quantity of 1.0 mass% addition 
at 24°C of diethylether is 0.1114 mL or 0.07953 g. 
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Fig. 1 Experimental apparatus. 


Table 1 Physical properties of Air & Diethylether[8,9] 


Work-fluid Air Diethylether C4H100 
Boiling point °C -193 34.8 
Density Liquid - 713.8 
kg/m? Vapor 1.181 3.308 
Heat of 

evaporation J/g : ice 

Cv J(gK) 0.7171 2.33:L 

à WK) 0.00257 0.129:L 


Experimental procedures 


The heat exchange cycle experiment is performed by 
the following procedures. 

(1) Exhaust the work-fluid which remains to test sec- 
tion by vacuum pump. 


(2) Air as a work-fluid of the temperature of 24°C and 
pressure 101.3 kPa (absolute humidity at this time is 
0.0112 + 0.0007 kg/kg*) is introduced into test section. 
Diethylether of specified quantity is thrown into test sec- 
tion by using digital burette. 

(3) A cooling surface and a heating surface are set to 
the experimental condition temperature. 

(4) The displacer piston starts to move to the upper 
side from the bottom side. The heat exchange cycle is 
continued. Then, the pressure variation reaches to the 
stable cycle. At this time, the pressure, heat transfer sur- 
face temperature, and work-fluid temperature of test sec- 
tion are measured, and it collects in a data logger. 


Experimental conditions 


Table 2 indicates the experimental conditions of this 
study. Temperature difference of heat source is set as 55- 
70K containing the boiling point of phase change materi- 
al. The high temperature heat source is set <80°C that 
aims as source temperature of hot spring. The low tem- 
perature heat source is maintained 10°C. The speed of a 
displacer piston is set as low-speed rotation N = 5-20 
r/min (ume = 0.2-0.9 m/s). The mixed rate to the work- 
fluid (air) of PCM (diethylether) is evaluated in A, = 0 ~ 
4 mass%. 

Fig.2 indicates the time history of displacer piston po- 
sition and velocity at this experimental condition. The 
position of a displacer piston is decided by the rotation 
angle regardless of number of revolutions. Sine wave 
operation is performed in the cycle process of a displacer 
piston, because a displacer piston on the actual engine is 
controlled by rotation angle. The displacer piston oper- 
ates by the computer controlled actuator. The heating 
process times are 1.5 seconds at 20 r/min, 3.0 seconds at 
10 r/min and 6.0 seconds at 5 r/min. 


Table 2 Experimental parameters and conditions 
45-60 step 5K 
55-70 step 5 K 


Temp. diff. bet. source AT K 
High temp. source T, °C 
Low temp. source T, °C 10 

Num. of rev. N r/min (Re) 5-20 (2784-8643) 


Mass composition ratio Aq mass% 0-4 


The cycle output is examined by using the cycle mean 
velocity on each number of rotations. The cycle mean 
velocity umc is calculated by the following equations. 

H = St xsin (1) 
_ N AHXWpp x Dpp 


u= (2) 
60 2xHxDpp 
1 0.57 
Unc =—f  udð (tp 20.1sec) (3) 
bor —0.57 
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Fig.2 Position and velocity of displacer piston on process. 


Fig.3 indicates the variation of the cycle mean velocity 
on the heat transfer surface with number of revolutions. 
The cycle mean velocity is related linearly to the number 
of rotations. 
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Fig. 3 Variation of cycle mean velocity with number of revo- 
lutions. 


The cycle mean velocity is estimated by the following 
equations as Reynolds number Re (using an averaged 


flow velocity of the process on each number of rotations). 


The characteristic length is set to 1/2 (L = 0.15 m) of heat 
transfer surface width, because the work-fluid flows into 
the center from both ends. The amount of additional 
phase change material is very little, so that the dynamic 
viscosity of air is used for calculation. 


_ Umel 


Re (4) 


Vair 
The quantity of heat which a work-fluid receives in 
heat exchange process is calculated by the following eq- 


uation with the application of a constant-volume change. 


0-6,C,[% T.-P) 6) 


In which, C,, and G, are the specific heat at constant 
volume and the mass of work-fluid. These values are 
calculated by using the thermophysical properties of 
composition. 

Cy = VairG 
(1- Add)x G 


+ Cv, Gy 
+ Add x Gł 


(6) 
(7) 


air 


Gy = air 


The quantity of heat in a cycle which a work-fluid 
receives is calculated by next integral equation. The ob- 
tained heat amount equals output per cycle. 


9 = fqao-w (8) 


The output improvement effect by the PCM addition is 
estimated by the output ratio. The output ratio is defined 
by next equation. 

w“ W gaa Waic (10) 
The temperature difference ratio is calculated by the 
following equation using PCM boiling point. 


T* {7 -T,)/(T,-T,) (11) 


Experimental Results and Discussions 


Indicator (AP-V ) diagram 

Fig.4 shows the relationship between the pressure dif- 
ference and the heating side volume (indicator AP-V dia- 
gram) at Aqa = 0 mass% (Air), N= 20 r/min, AT = 50K. 
The cycle pressure changes to the high pressure from the 
start of experiment. The stable cycle focusing on OkPa is 
observed at n= 22. 
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Fig. 4 Relationship between the differential-pressure and the 
heating side volume (44a = 0 mass%, N = 20 r/min, 
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Fig. 5 Variation of the output per cycle with the times of 


driving cycle. (aa = 0 mass%, AT = 45K) 


Variation of the output per cycle with the number of 
cycle is shown in Fig.5. The output per cycle is calcu- 
lated by cycle integration of each cycle. The experimen- 


tal data are indicated by the plots in this figure. The 
curved lines show the approximation values that are ob- 
tained from experimental data. Output reaches the stable 
value by driving for more than 80 seconds in spite of the 
condition. Therefore the data of drive time after 120 
seconds as the stable cycle are used in this study. 


Cycle output by using air as work-fluid 


Fig.6 shows the variation of the output per cycle with 
the number of revolutions at Aaa = 0 mass% (air). Cycle 
output decreases with the increasing of the number of 
rotations. Because the decreasing of heat exchange time 
is brought by the increasing of number of rotations. 

Fig.7 shows the variation of the output per cycle with 
the temperature difference of heat source at Ag = 0 
mass% (air). Cycle output increases with the increase of 
temperature difference of heat source. 

The experimental data of cycle output by using air 
(Aaa = 0 mass%) as work-fluid is correlated as following 
equations. The values obtained by the equations are equal 
to the experimental data within + 10 %. 


=a 0 (5<N20) (12) 


W, 


air 


a =10.1AT — 230 (45 <AT<60) (13) 
b=0.0405AT? —3.81AT +97.0 (45SATS60) (14) 
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Fig. 6 Variation of the output per cycle with the number of 


revolutions. (A44 = 0 mass%) 
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Fig. 7 Variation of the output per cycle with the temperature 


difference of heat source. (Aaa = 0 mass%) 


Cycle output by using air-PCM mixture as work-fluid 


Fig.8 indicates the relationship between the output per 
cycle and the cycle time or the number of rotations at Aqa 
= 0.5mass%. Output per cycle is decreased with the in- 
creasing of number of rotations. 

Fig.9 indicates the variation of output ratio with tem- 
perature difference ratio at N=15 r/min (Re = 6870). The 
plots in the figure indicate experimental value. The lines 
in the figure indicate the value calculated by the follow- 
ing experimental correlation equation. Output ratio shows 
the high value at the range of the temperature difference 
ratio AT =0.8-1.2. 
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Fig. 8 Relationship between the output per cycle and the cycle 
time or the number of rotations. (444 = 0.5 mass%) 
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Fig. 9 Relationship between the output ratio and the heat 
source temperature difference ratio. (N = 15 r/min) 
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Fig. 10 Variation of the output ratio with Reynolds number or 
the number of rotations. (AT = 55K) 


Fig.10 shows the relationship between the output ratio 
and the Reynolds number or the number of rotations at T“ 
= 1.2. It is observed that the high output ratio range is Re 
= 4945-6870 (N = 10-15 r/min). 

Fig.11 indicates the relationship between the output 
ratio and PCM additional ratio at Re = 6870 (N = 15 
r/min). The high output ratio is indicated at the range of 
Agq= 0.5-1 mass%. 
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Fig. 11 Relationship between the output ratio and the PCM 
additional ratio. (N = 15 r/min) 


Output ratio W is correlated by Ag, A T * and Re. 
Therefore, the following equation is obtained. Applicable 
range of the equation is Re = 4945~8643 (N = 10~20 
r/min), AT’ = 0.8~1.2 (AT = 45~55K) and Aa = 0.5~3 
mass%. 
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Conclusions 


Experimental study of the heat exchange cycle by us- 
ing the work-fluid including PCM (diethylether) is car- 
ried out, and the following conclusions are obtained. 

Output improvement by PCM addition work-fluid was 
observed in the heat exchange cycle. 

The optimum range of condition (mass composition 
ratio of PCM, speed and temperature) for the cycle out- 
put improvement was indicated . 

The useful experimental correlation equation was 
suggested by the cycle output ratio. 
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